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Abstract 

The spherical infall model first developed by Lemaitre and Tolman was modified 
in order to include the effects of a dark energy term. The resulting velocity-distance 
relation was evaluated numerically. This equation, when fitted to actual data, per- 
mits the simultaneous evaluation of the central mass and of the Hubble parameter. 
Application of this relation to the Local Group, when the dark energy is mod- 
eled by a cosmological constant, yields a total mass for the M31-Milky Way pair 
of (2.5 ± 0.7) X 10^^ Mq, a Hubble parameter Hq = 74 ± Akms'^Mpc-^ and a 
1-D velocity dispersion for the flow of about 39 kms^^. The zero- velocity and the 
marginally bound surfaces of the Local Group are at about 1.0 and 2.3 Mpc re- 
spectively from the center of mass. A similar analysis for the Virgo cluster yields a 
mass of (1.10 ±0.12) x 10^^ M© and Hq = 65 ± 9 fcms-^Mpc-^ The zero-velocity is 
located at a distance of 8.6±0.8 Mpc from the center of the cluster. The predicted 
peculiar velocity of the Local Group towards Virgo is about 190/cms~^, in agree- 
ment with other estimates. Slightly lower masses are derived if the dark energy is 
represented by a fluid with an equation of state P = we with w = —2/3. 
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1 Introduction 



New and high quality data on galaxies belonging to nearby groups have im- 
proved considerably estimates of masses and mass-to-light ratios (M/L) of 
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these systems. Searches on the POSS II and ESO/SERC plates (Karachent- 
seva & Karachentsev 1998; 2000) as well as "blind" HI surveys (Kilborn et 
al. 2002) lead to the discovery of new dwarf galaxies, increasing substantially 
their known population in the local universe. Moreover, in the past years, us- 
ing HST observations, distances to individual members of nearby groups have 
been derived from magnitudes of the tip of the red giant branch by Karachent- 
sev and collaborators (Karachentsev 2005 and references therein), which have 
permitted a better membership assignment and a more trustful dynamical 
analysis. 

Previous estimates of M/L ratios for nearby groups were around 170Mq/ Lb.q 
(Huchra & Geller 1982). However, virial masses derived from the aforemen- 
tioned data are significantly smaller, yielding M/L ratios around 10-30 Mq/ Lb,q- 
If these values are correct, the local matter density derived from nearby groups 
would be only a fraction of the global matter density (Karachentsev 2005). 
However, for several groups the crossing time is comparable or even greater 
than the Hubble time and another approach is necessary to evaluate their 
masses, since dynamical equilibrium is not yet attained in these cases. 

Lynden-Bell (1981) and Sandagc (1986) proposed an alternative method to 
the virial relation in order to estimate the mass of the Local Group, which can 
be extended to other systems dominated either by one or a pair of galaxies. 
Their analysis is essentially based on the spherical infall model. If the motion 
of bound satellites is supposed to be radial, the resulting parametric equations 
describe a cycloid. Initially, the radius of a given shell embedding a total mass 
M expands, attains a maximum value and then collapses. At maximum, when 
the turnaround radius Rq is reached, the radial velocity with respect to the 
center of mass is zero. For a given group, if the velocity field close to the main 
body, probed by satellites, allows the determination of Rq, then the mass can 
be calculated straightforwardly from the relation 

M = ^ (1) 



where To is the age of the universe and G is the gravitational constant. 

Data on the angular power spectrum of temperature fluctuations of the cosmic 
microwave background radiation derived by WMAP (Spergel ct al. 2003) and 
on the luminosity-distance of type la supernovae (Riess et al. 1998; Perlmutter 
et al. 1999), lead to the so-called "concordant" model, e.g., a fiat cosmolog- 
ical model in which = 0.3 and Vty — 0.7. The later density parameter 
corresponds to the present contribution of a cosmological constant term or a 
fluid with negative pressure, dubbed "quintessence" or dark energy. The ra- 
dial motion leading to the aforementioned M = M{Rq,Tq) relation neglects 
the effect of such a term, which acts as a "repulsive" force. This repulsive 
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force is proportional to the distance and its effect can be neglected if the zero- 
velocity surface is close to the center of mass. Turnaround radii of groups are 
typically of the order of 1 Mpc (Karachentsev 2005), while the characteris- 
tic radius at which gravitation is comparable to the repulsion force is R* = 
I.IM}^^ Mpc, where M12 = M/(1O^^M0) and the Hubble parameter Hq was 
taken equal to 70 kms~^ Mpc~^ . This simple argument suggests that the effect 
of the cosmological term can not be neglected when deriving masses from the 
M = M{Ro,To) relationship. 

In this paper we revisit the velocity-distance relationship when a dark en- 
ergy term is included in the dynamical equations and calculate the resulting 
M = M{Rq, Hq) relation. The presence of the dark energy term has been 
also invoked as a possible explanation for the smoothness of the local Hubble 
flow (Chernin 2001; Teerikorpi, Chernin and Baryshev 2005) being a further 
reason to investigate its effects on the M — M(Rq,Tq) relation. For a cosmo- 
logical density parameter associated to the "vacuum" energy Qy — 0.7 (our 
preferred solution), numerical computations indicate that the "zero-energy" 
surface, beyond which galaxies will never collapse onto the core, is located at 
about 2.3it!o- In order to illustrate our results, some apphcations are made to 
the Local Group and the Virgo cluster. As we shall see, values of the Hub- 
ble parameter resulting from fits of the actual data to the velocity-distance 
relation including a cosmological term, are in better agreement with recent 
estimates then those derived from the relation obtained either by Lyndell-Bell 
(1981) or Sandage (1986). In Section 2 the relevant equations are introduced, 
in Sections 3 and 4 the results are applied to the Local Group and the Virgo 
cluster and finally, in Section 5 the concluding remarks are given. 



2 The velocity-distance relation 

The evolution of a self-gravitating zero-pressure fluid with spherical symmetry 
was flrst considered by Lemaitre (1933) and Tolman (1934). The Lemaitre- 
Tolman model describes quite well the dynamics of an extended halo around 
a bound central core, asymptotically approaching a homogeneous Priedmann 
background. In this situation, three main distinct regions can be distinguished: 
i) the central core, in which the shell crossing has already occurred, leading to 
energy exchanges which transform radial into transverse motion; ii) the zero- 
velocity surface, boundary which separates infalling and expanding bound 
shells and iii) the "marginally" bound surface (zero total energy), segregat- 
ing bound and unbound shells. Density proflles resulting from the Lemaitre- 
Tolman model were examined by Olson & Silk (1979) and application of this 
model to the velocity held close to the Virgo cluster were made by Hoffman 
et al. (1980), TuUy & Shaya (1984), Teerikorpi et al. (1992) among others. 
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If displacements of galaxies, here associated to the outer halo shells, develop 
mainly at low redhsifts when the formation of the mass concentration around 
the core is nearly complete (see, for instance, Peebles 1990), then the equation 
of motion for a spherical shell of mass m, moving radially in the gravitational 
field created by a mass M inside a shell of radius R, including the dark energy 
term is 

d'^R GM (1 + 3W) ^ rr2D/«0x3(l+«,) 



where M — An J^r'^pmdr and a is the scale parameter (the present value is 
taken as oq = 1). The latter satisfies the Hubble equation 

(^)^ = i^o[^^.(^)^ + ^^.(^)^(^+-)] (3) 



Here, the common assumption that the dark energy can be modeled as being 
a fluid with an equation of state P = we was adopted and in the two equations 
above, the dependence of the dark energy on the scale parameter was obtained 
by solving the energy conservation for such a component. Eq. (2) is intended 
to describe the motion of shells in the halo, excluding the central region where 
shell crossing effects have probably already occurred. 

Defining the dimensionless variables y = R/Rq, t = tHo and x = a/ao, eqs. 
(2) and (3) can be rewritten as 

|| = -il| + B..-3(-.l (4) 



and 



dx /^m /r\ 



where we have introduced the parameters A = 2GM/{HqRq) and B = (1 + 
3w)Q^. These equations were solved numerically by adopting the following 
procedure. For a given rcdshift, the initial value of the scale parameter is 
derived as well as the corresponding instant of time from the Hubble equation. 
If initially, at high redshifts (here taken around z ~ 100), the dark energy term 
is negligible, then using a Taylor expansion of the standard Lemaitre-Tolman 
solution, when the angle parameter 9 << 1 (sec, for instance, Peebles 1980), 
the initial values of y and its derivative dy/dr can be estimated. For a given 
value of w, the parameter A is varied until the condition defining the zero- 
velocity surface, e.g., dy /dr — at y — 1 is satisfied. For the case ty = — 2/3, 
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we have obtained A = 3.414 and for the particular case w = —1, representing a 
cosmological constant, A — 3.658. Therefore, the mass inside the zero- velocity 
radius Rq is 

M = 1.827^^^ = 4.1 X lO'^h^Rl Mq (6) 



where h = Ho/[100 kms~^Mpc~^) and Rq is in Mpc. In the case w = —2/3, 
the numerical coefficient is slightly smaller (1.705 instead of 1.827). Comparing 
with eq. (1), we notice that the inclusion of the dark energy term represents, 
for a given Rq, an increase of about 28-38% on the mass derived by such a 
procedure. 

Once the parameter A is known, the velocity-distance relation, v = v(R), for 
different shells at a given time is obtained by varying their energy. Shells with 
negative energy will expand, halt and fall back toward the center, while shells 
with positive energy expand forever, according to the aforementioned char- 
acterization of regions ii and Hi. At a given time, there is a critical energy 
Ec which defines the zero-velocity radius. Shells having E < Ec have already 
crossed the turnaround point and are collapsing. Consequently, they have neg- 
ative velocities. Shells with E > E^ are still expanding and thus have positive 
velocities. 

For the case w = — 1, which gives a good representation of actual data, as 
we shall see in the next sections, the resulting numerical values are quite well 
fitted by the relation 



with n = 0.865. Notice that from the condition v{Rq) = 0, eq. (6) is recovered. 
It is worth mentioning that solutions with negative energies are possible up to 
it! = 2.30it!o, which defines the marginally bound surface. Beyond this critical 
radius, only unbound shells exist. The solution for w = —2/3 differs only 
marginally and will not be considered in further analyses. 

For comparison, the velocity- distance relationship derived for the case ^2^ = 
is 

v{R) = -1.038(^)^/' + 1.196//oi? (8) 



which has a slightly fiatter dependence on the distance than the precedent 
equation. 
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3 The Local Group 



An immediate application of the velocity-distance relationship is the determi- 
nation of the Local Group mass, concentrated mainly on M31 and the Milky 
Way, as well as the Hubble parameter itself as we shall see later. Moreover, 
the velocity-distance relation gives also an indication of the dispersion of the 
peculiar velocities over the Hubble flow. As we have already mentioned, the 
local velocity dispersion is known to be quite small (Giraud 1986; Schlegel et 
al. 1994), a fact referred usually as the "coldness" of the local flow. An in- 
vestigation of the dynamics of the Local Group and its environment by using 
numerical simulations was performed by Governato et al. (1997), who con- 
cluded that cold dark matter models (f)^ = lor flm = 0.3) are unable to 
produce candidates embedded in regions having "cold" flows. 

As a flrst step, we have searched to check the accuracy of eqs. (6) and (7) by 
using numerical simulations. We have performed N-body simulations using the 
adaptive particle-particle/particle- mesh (AP'^M) code HYDRA (Couchman 
et al. 1995), with cosmological parameters h = 0.65, = 0.3, = 0.7 
and (Tg = 0.9. The simulation was performed in a periodic box of side 30h~^ 
Mpc including 256^ particles, corresponding to a mass resolution of 2.05 x 
10^ Mq. The simulation started at z = 49 and ended at the present time. Halos 
were initially detected by using a friends-of-friends (FOF) algorithm and, in a 
second step, unbound particles were removed by an iterative procedure. Thus, 
all halos in our catalog are gravitationally bound objects. For further details, 
the reader is referred to Peirani et al. (2004). 

In our mock halo catalog, several examples of pairs with physical characteris- 
tics similar to the MW-M31 pair can be found. However, most of them have 
nearby (within 3-5 Mpc) halos of comparable mass or even higher, which per- 
turb considerably the velocity field. Here, for illustration purposes only, we 
consider one case in which the main halos have masses respectively equal to 
9.85 X 10^"*^ and 6.48 x 10^"*^ solar masses and are separated by a distance of 
0.69 Mpc, parameters comparable to those of the M31-MW pair. About 73 
subhalos were detected within 3 Mpc of this pair, but satellites situated at 
distances larger than 2 Mpc are clearly perturbed by other nearby structures, 
since the velocity dispersion increases considerably at those distances. Veloc- 
ities and distances with respect to the center of mass were computed for all 
these objects. 

A simple fitting of these simulated data gives Rq = 0.98 ± 0.20 Mpc and 
using eq. (6) with h = 0.65 (adopted in the simulations) a total mass of 
(1.60 ± 0.32) X 10^^ Mq is obtained. In spite of the excellent agreement with 
the actual total mass of the pair, this result is somewhat fortuitous since the 
velocity dispersion of the simulated satellites is rather high. In practice, besides 
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Fig. 1. Simulated velocity-distance data (diamonds) and best fit to the v=v(R) 
relation (solid curve), corresponding to M = 1.48 x 10^^ Mq and a = 73 km/s. The 
Hubble parameter was held constant, h = 0.65. 

random motions, errors in distances or velocities increase considerably the 
uncertainty in the determination of the zero- velocity radius. These difficulties 
can be alleviated by searching the best fit of the v — v(R) relation to data. 
By varying the mass in order to minimize the velocity dispersion and giving a 
higher weight to the inner satellites, one obtains M = (1.48±0.30) x 10^^ Mq. 
This result seems to be more confident to evaluate the uncertainties of the 
method. It is worth mentioning that quoted errors are estimates based on the 
spread of values derived from the fitting procedure and not formal statistical 
errors. In figure 1 we show the simulated velocity-distance data for satellites 
with R < 1.8 Mpc and the best fit solution for cq. (7). The derived 1-D 
velocity dispersion for this simulated data is 73 kms^^, a value lower than 
that derived by Governato et al. (1997) but in quite good agreement with the 
ACDM simulations by Maccio et al. (2005), who obtained a velocity dispersion 
of about 80 kms~^ within a sphere of 3 Mpc radius (see their Fig. 3). A fiat, 
A dominated cosmology is able to produce fiows on scales of few Mpc around 
field galaxies "colder" than pure dark matter models, but somewhat higher 
than values derived from actual data, as we shall see below. 

Recent data on neighboring galaxies of the Local Group were summarized by 
Karachentscv ct al. (2002), who have estimated Rq = 0.94 ± 0.10 and derived 
from eq. (1) a total mass of 1.3 x 10^^ ^-j^g M31/MW pair. 

Here, eq. (7) was fitted to the data by Karachentsev et al. (2002), but varying 
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Fig. 2. Velocity and distance data (circles) for satellites of the M31-MW system 
(Karachentsev et al. 2002) and the best fit to the v=v(R) relation, corresponding 
to M = 2.5 X 10^2 ^ and h = 0.74 



now both the mass and the Hubble parameter in order to minimize the velocity 
dispersion. We have obtained h = 0.74 ±0.04 and M = (2.5 ±0.7) x 10^^ M©, 
where the quoted errors are again estimates based on the uncertainties of 
the fitted parameters. Figure 2 shows data points and the velocity-distance 
relation defined by the previous parameters. Had we used eq. (8) instead of 
eq. (7) in the fitting procedure, a similar result for the mass would have been 
obtained, but with a /li^/ier Hubble parameter, e.g., h — 0.87 ±0.05. We shall 
return to this point latter. 

The zero- velocity radius is located at 1.0 ±0.1 Mpc, which is about 10% higher 
than the value estimated by Karachentsev et al. (2002) and the "marginally" 
bound surface (zero-energy) is at a distance of about 2.3 Mpc. Prom our fit 
it results a 1-D velocity dispersion of 39 kms~^, which should be compared 
with the value of 73 km/s found from our simulated data for a similar scale. 
Maccio et al. (2005) have also revisited the "coldness" of the Hubble fiow and, 
according to their final results, within a sphere of radius 3 Mpc, the expected 
velocity dispersion is 38 kms~^, in rather good agreement with our figure. 
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Fig. 3. Velocity and distance data for galaxies with Virgocentric distances in 
the range 3.5 < R < 15 Mpc and the best fit to the v=v(R) relation, for 
M = 1.1 X 10^5 and h = 0.65 

4 The Virgo Cluster 

Dynamical models for the Virgo cluster based on the Lemaitre-Tolman model 
were, for instance, developed by Hoffman et al. (1980). They have modeled the 
projected velocity dispersion as a function of the angular distance and, from 
comparison with data, derived a mass of (4.0 ± 1.0) x 10^'^h~^ Mq contained 
inside a sphere of 6° radius, which corresponds approximately to the central 
relaxed core of the cluster. Using the virial relation, TuUy & Shaya (1984) 
obtained a mass of (7.5 ± 1.5) x 10^^ Mq for this central core. More recently, 
Fouque et al. (2001) using the Lemaitre-Tolman model derived a mass of 
1.3 X 10^^ Mq inside a radius of 8°. 

In this section, the derived velocity-distance relation (eq. 7), including effects 
of the cosmological constant, will be applied to galaxies outside the inner core 
of the Virgo cluster. 

Galaxies with Virgocentric distances higher than 1.7 Mpc, corresponding ap- 
proximately to the core radius, and less than 15 Mpc, were selected from 
the list by Teerikorpi et al. (1992), constituting a sub-sample of 27 objects. 
Distances derived from the Tully-Fisher relation were taken from the afore- 
mentioned source. Heliocentric velocities were corrected with respect to the 
Local Group according to the prescription by Courteau and van den Bergh 
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(1999). Velocities and distances with respect to Virgo center were calculated 
by assuming a distance to the cluster of 16.8 Mpc (Tully & Shaya 1984) and 
an observed velocity of 967 km/s (Kraan-Korteweg 1981). 

The best fit of eq. (7) to data gives h = 0.65 ± 0.09 and M = (1.10 ± 0.12) x 

10^^ Mq, corresponding to a velocity dispersion (Tid = 335 km/s for the flow. 
Figure 3 shows the velocity-distance data for the galaxies of our sample and 
the theoretical v = v(R) relation computed with the derived parameters. The 
higher mass derived in the present analysis confirms some early results based 
on models of the velocity field in the vicinity of the Virgo cluster using the 
Lemaitre-Tolman equations, as those performed by Tully & Shaya (1998) and 
Fouque et al. (2001). 

The resulting zero-velocity surface is located at Rq — 8.6 ± 0.8 Mpc, which 
at the assumed distance corresponds to an angle 6zv = 29°. This should be 
compared with the analyses by Hoffman et al. (1980), Tully & Shaya (1984) or 
Teerikorpi et al. (1992), who estimated 9zv = 27", 28" and 25.8" respectively. 

Eq. (7), with the above parameters, predicts that at the level of the Local 
Group the observed velocity should be 988 km/s, which compares quite well 
with the observed value. This indicates that the projection of the peculiar 
velocity of the Local Group in the Virgo direction is about 190 km/s, compat- 
ible with the values found by de Freitas Pacheco (1985), Tamman & Sandage 
(1985) and Federspiel et al. (1998). The marginally bound surface is located 
at about 19.8 Mpc, implying that the Local Group is bound to the cluster. 



5 Conclusions 

The contribution of a dark energy term in the mass-energy budget of the 
universe seems to be well established at the present time. In this study, the 
usual velocity-distance relation based on the Lemaitre-Tolman model, was 
revisited in order to include effects due to such a cosmological term. 

The dynamical equations were solved numerically and the relation M — 
M(Ho, Rq), defining the mass inside the zero- velocity surface was recalculated. 
For a given Rq, the resulting masses are about 35-38% higher with respect to 
the original relation derived from the Lemaitre-Bondi model {fly = 0), if the 
dark energy is modeled by a fluid with an equation of state P/e — w — —2/3 
or by a cosmological constant respectively. 

The resulting v = v(R) relation {w = —1 case) was applied to the Local Group 
and to the Virgo cluster. From the best fitting procedure, a mass of (2.5±0.7) x 
10^^ Mq was derived for the M31-MW system, which corresponds to a M/L 
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ratio of about 25 Mq/Lb,q- The zero- velocity surface and the "marginally" 
bound surface are located at about 1.0 ±0.1 Mpc and 2.3 Mpc respectively. 
The later defines in fact the boundaries of the Local Group. 

More sophisticated analyses of the Local Group based on the action principle 
were performed by Peebles (1990). In this method, orbits are reconstructed 
in order to reproduce the observed radial velocities, maintaining the action 
stationary. Peebles (1990) preferred solution implies a mass of 6.4xlO^^M0 
for the Local Group, a factor 2.5 higher than our result, but a comparison of 
these results could be meaningless since a different cosmology {VL = 0.1) was 
adopted in that work. 

The same procedure apphed to the Virgo cluster gives a mass of about (1.10 ± 
0.12) X 10^^ Mq inside a radius of 8.6 ±0.8 Mpc. This mass is higher than the 
virial value (Tully & Shaya 1984) or values derived from velocity dispersion 
profiles computed from models based on the Lemaitre-Tolman model (Hoff- 
man et al. 1980). However, recent studies of the velocity field based on the 
Lemaitre-Tolman formulation also lead to masses compatible with our result. 
Our velocity-distance relation predicts a peculiar velocity of the Local Group 
towards the Virgo cluster of about 190 km/s, compatible with different esti- 
mates. The "marginally" bound surface encloses the Local Group, which in 
the future may attain its zero-velocity surface and then fall onto the cluster. 
However this prediction is based on a model in which the Local Group is con- 
sidered as an "isolated" system. In reality the motion of the Local Group is 
much more complex, being strongly affected by large mass concentrations in 
the direction of Hydra- Centaurus (Great Attractor) and the Shapley super- 
cluster. Some numerical simulations predict that within ~ 30 Gyr the Local 
Group will get closer to the Virgo center (in comoving coordinates) , but then 
will be pulled away(in physical coordinates) due to the accelerated expansion 
of the Universe (Nagamine & Loeb 2003). 

The introduction of a cosmological constant term modifies the velocity-distance 
relation in comparison with that derived from the Lemaitre-Tolman model. 
Nevertheless both descriptions of the velocity field near the Local Group or 
the Virgo cluster yield masses comparable to within a factor of two. This is 
probably due to the fact that errors still present in distance estimates mask 
differences between both models. However, when searching for a best fit of 
both models to data, there is a substantial difference in the resulting Hubble 
parameter. The Lemaitre-Tolman model requires h in the range 0.87-0.92 in 
order to fit adequately the Local Group and Virgo data respectively, whereas 
eq.(7) requires h in the range 0.65-0.74, more consistent with recent determi- 
nations and with the "concordant" model. In this sense, the inclusion of the 
cosmological constant in the v(R) relation seems to improve the representation 
of actual data. 
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